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Abstract 

The problem of the transformation of microscopic information to the macro- 
scopic level is an intriguing challenge in computational neuroscience, but also 
of general mathematical importance. Here, a phenomenological mathematical 
model is introduced that simulates the internal information processing of brain 
compartments. Synaptic potentials are integrated over small number of realis- 
tically coupled neurons to obtain macroscopic quantities. The striatal complex, 
an important part of the basal ganglia circuit in the brain for regulating motor 
activity, has been investigated as an example for the validation of the model. 



1. Introduction 

The brain nuclei, as parts of complex brain networks, are comprised by 
different types of inter- and projection neurons within subnetworks of highly 
complex structure. Beside the anatomical formation, functional properties in 
integration and processing of neural information become adapted during the 
development of these regions until adolescence. In this study, adult brain nuclei 
are considered after their final development, as part of information processing 
pathways. 

Experiments on the topology of brain regions suggest that these structures 
can be categorized as heterogeneous media. There are various mathematical 
methods describing dynamical processes in heterogeneous media, such as asymp- 
totic analysis, homogenization, and integral equations (Pavliotis and Stuart |17j . 
Haken [T5], E and Engquist [8j, Levin and Chaop^). There are traditional net- 
work approaches on the effects of properties such as the local connectivity of 
neurons on the striatal function (Wickens et al. [SI])- These studies concen- 
trate merely on the network structure of a particular brain region e.g. striatum 
and even then do not contain some important features such as the role of large 
spiny neurons in the striatal function. Large spiny neurons in the striatum are 
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cholinergic and are of great importance for processes involving addiction and 
food intake (Rada et al. [T^, Avena et al. [T]). 

In the present study, simplified but essential physiological processes have 
been chosen as the biological foundation for the mathematical model. The 
reader is advised to consider the appendix for a brief introduction into the 
biological terminology. 

The course of information processing in brain regions depends on the propa- 
gation of synaptic potentials - which are averaged postsynaptic potentials - along 
the subnetworks of these regions. Thus, understanding the physiology of this 
integration process requires the comprehension of the neuronal architecture of 
the brain region of interest. The regional subnetworks consist of different types 
of neurons and characteristic connections among them that can be obtained 
by experiments. In this study the intraneuronal connections are subdivided 
into the global, and ultrastructural morphology. Golgi and immunohistochem- 
ical methods provide optimal frameworks for obtaining the values of the two 
morphological dimensions experimentally. The immunohistochemical studies 
reveal the appearance of certain neurotransmitter systems in a structure; Golgi 
studies provide information about the form and type of neurons expressing the 
obtained neurotransmitter systems and on the morphological dimensions of the 
brain region in terms of the characteristic connections between different types 
of neurons in a brain region. In general, the structural topology of the brain 
region is then obtained. The knowledge about the appearance of different types 
of neurons and the intraneuronal synaptic interactions inside a nucleus allow us 
to abstract its structure through networks. 

One way of representing such networks is through their realization in con- 
tinuous spaces. The attendance of a high amount of neurons supports the idea 
of the appropriateness of such an approach for this purpose. 

Because of the importance of the global and ultrastructural morphology for 
the information integration process in a nucleus and the dependency of their 
values on the spatial distribution of neurons, they shall be modelled as spatial 
variables parametrized by the neuron (neurotransmitter) type. 

To provide a general framework for the regional neuronal activity, we intro- 
duce an integral operator averaging the information on two refinement scales. 
The idea is to cover the space including the nucleus with copies of discrete 
fundamental domains consisting of neuronal assemblies which we call n-cells. 
Roughly spoken, n-cell patches are finite networks of neurons. The ultrastruc- 
tural morphology of a nucleus in terms of the synaptic connectivities between 
different neurons provides information on the local connectivities in the n-cells, 
which are embedded as edges of neuronal network inside the n-cells. Within the 
information on the adherence functions of n-cells (global morphology variables) 
and the distribution of neurons, the averaging of synaptic potentials across the 
network of nuclei is completed. The main advantage of this method is in its 
appropriateness for applications on various brain regions. 

Assembling a suitable number of n-cells by considering the large-scale inter- 
action between the n-cells (adherence of n-cells), we obtain information on the 
combined influence of several synaptic inputs on macroscopic neuronal activity 
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Figure 1: Composition of n-cells 



(Fig. [T|). One trivial example is the subthalamic nucleus which consists of only 
one class of interneurons. In this case the number of n-cells under consideration 
is one (Parent and Parent [16 ). 

The aim of this study is to present a mathematical technique for averaging 
microscopic information along complex network structures such as the networks 
of single brain regions. The technique is applied on the corpus striatum - a part 
of the basal ganglia circuit which plays an important role in the regulation of 
motor behaviour- and shown to mimic faithfully the qualitative oscillatory be- 
haviour of this brain region. The general formulation of the averaging operator 
suggest its suitability for other discrete multiscale problems, especially in the 
research area of material sciences. 

2. Model Statement 

Statement:. The integration of the electric activity in a brain nucleus (infor- 
mation processing) depends on the distribution of its comprising neurons, the 
ultrastructural morphology as a local variable quantity and the global quantity 
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in terms of the composition of neuronal fundamental domains (n-cells), that 
describes the topology of the brain region. 

The concept of the n-cells is very essential for the understanding of the 
averaging transformation. The idea is to path up a region with multiple copies 
of discrete fundamental domains called n-cells. The n-cells compose the whole 
brain region of interest by applying the global morphological properties and 
contain the local morphological properties as finite graphs. Hence, they unify 
the morphological dimensions into one notion. 

Definition 1: A n-cell is a neuronal assembly represented by a finite graph 
characterized by the local ultrastructural morphology. Within proper composi- 
tion of the n-cells, as discrete fundamental domains, they patch up the region 
of interest. 

Let fl C a compact, path-connected set, be the brain region of interest. 

Let C ri, i G / be a skeleton representation of the finite graph of a n-cell, 
and / be the finite index set with #/ = #{n-cells}. Furthermore, let x € fij be 
the variable describing the network position of a neuron relative to the n-cell 
fli] and y the continuity variable describing the spatial position of a neuron 
embedded as points in as a subset of the complete space R"^. {z} denotes the 
finite set of neurotransmitter families that appear in the brain region of interest. 
By allocating any neuron to a neurotransmitter, z parametrizes the proposed 
averaging. 

Definition 2: A synaptic potential u^{t) is the sum of all Hodgkin-Huxley 
postsynaptic potentials of a neuron of type z at position x inside a n-cell. 

yi€ I , Wz , Wxefli : < : [0,T] ^ R . 

To average the synaptic potentials u^{t) in n-cells, the local connectivity 
function : fli ^ R+ has to be characterized. By introduction of the global 
connectivity function : — > R+, describing the composition of n-cells, we will 
cover n will n-cell patches. We write x (Xi) for the |/|-vector of composition 
mappings. The normalized distribution function : f2 — > [0, 1] as the kernel of 
the averaging operator will complete the information required for the integration 
of synaptic potentials across the network of brain regions. These functions 
(■0^) Xi) P^) are obtained from Golgi- and immunohistochemical studies. 

The continuous realization of the network structure of the brain nuclei and 
the idea of transforming the microscopic, synaptic information to the macro- 
scopic level compatible with the internal structure of a brain region suggest the 
construction of an averaging integral operator along networks. 

Definition 3: A compartment is a multiple C := (fi, p^, x, V'^; -^) which 
characterizes a brain nucleus by its ultrastructural and global morphology. 
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Definition 4: Let C be a compartment, then the averaged potential v{t) of 
synaptic potentials u^{t) across the network of a single brain nucleus is repre- 
sented by: 



g : RI^I R+ is the proper averaging parameter estimated by immunohisto- 
chemical experiments. It provides a combinatorial dimension for the composi- 
tion of the n-ccUs in a compartment. 

The universality of the concept of a compartment provides the possibility 
of applying the averaging transformation of this type for a wide category of 
network multiscale problems. 

3. Numerical Simulations 

The aim of this section is the numerical investigation of the dynamical be- 
haviour of a compartment by known synaptic dynamics (in terms of electric 
activity u^(t)) of single neurons. The corpus striatum is used as an non-trivial 
example to illustrate the efficiency of the averaging operator method. The com- 
plexity of the striatal structure on the one hand, and the importance of this brain 
region as a substrate of the basal ganglia in the regulation of motor activity and 
neurological diseases on the other hand, make it a proper representative for such 
investigations. First, the main morphological and ultrastructural properties of 
this compartment are summarized which are required to comprise the integral 
function. These properties include the spatial distribution of the different neu- 
ral populations (p^), the neurochemical classification (z), and the intraneuronal 
connections ((V'^j Xij 5))- Then, the simulation results are discussed and com- 
pared with the experiments. 

A category of proper experiments for the validation of the numerical results 
is represented by the local field potential (LFP) studies. Thereby, a signal is 
recorded using a low impedance extracellular microelectrode, placed sufficiently 
far from individual local neurons to prevent any particular cell from dominating 
the electrophysiological signal. This signal is then low-pass filtered, cut off at 
« SOOiJz, to obtain the local field potential (LFP). The low impedance and 
positioning of the electrode allows the activity of a large number of neurons to 
contribute to the signal. The unfiltered signal refiects the sum of action poten- 
tials from cells within approximately 50 — 350/iTO from the tip of the electrode 
(Legatt et al. [H!). 

The frequency of the LFP signals in corpus striatum will then be compared 
with those appeared as the results of the numerical simulation of the averaging 
transformation by given morphological properties. 
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3.1. The Corpus Striatum 

The neuronal populations of the striatum could be divided into four classes: 
Spiny projection neurons (about 96% of the whole neural population) are GABAer- 
gic neurons, which get external inputs from cortical areas, thalamus and sub- 
stantia nigra pars compacta (Smith et al. [2(1). These neurons also get inputs 
from the dopaminergic, GABAergic and acetylcholinergic interneurons. Large 
aspiny neurons (cholinergic), GABAergic interneurons and dopaminergic neu- 
rons comprise the rest population of striatal neurons. The cholinergic interneu- 
rons get external inputs from thalamus and substantia nigra pars compacta; 
and internal inputs from the GABAergic interneurons. They project to the 
dopaminergic and projection neurons. The action of the cholinergic interneu- 
rons is assumed to have opposite effects on the action of the dopaminergic 
neurons on the whole network. The cholinergic neurons inhibit the activity of 
the projection neurons of the direct pathway (mostly Di system) and disinhibit 
the activity of the indirect pathway neurons (Z?2 system) (DiFiglia [Tj , Graybiel 
et al. [12] , DiFiglia et al. [5] , DiFigHa and Carey [6] , Gerfen [10] , Betarbet et al. 
[5], Flores-Hernandez et al. [S], Graveland et al. [TT], Bennett and Wilson [5]). 
A schematic representation of the microcircuitry of corpus striatum is shown in 
Fig.|2] 

MATLAB has been used to simulate the averaging across the network of 
striatum. An overall number of 6400 neurons has been used for the simulation 
of averaged potentials in striatum. The density distribution function p is repre- 
sented in Figure 3. The local and global connectivity function ^ and x defined 
the necessary computational parameters based on Figure 2. We have induced 
a single excitatory signal from cortical projection neurons to the cholinergic 
neurons in striatum to simulate the simplest integration process. 

It reveals that the activity transmission induced by a constantly activated 
neurons is approximately radial and oscillatory. The activity of a single neuron, 
forced by a constant input, produces activation- waves on the striatal populations 
(Fig. §. 

Such neurons are simplified versions of the tonically active cholinergic in- 
terneurons. Activations of neurons of other classes suggest similar activation 
waves across the striatal populations. We observe that the simulated oscil- 
lations of v{t) (« bQHz) are qualitatively correlated with the LFP-studies of 
basal ganglia which also suggest oscillations in control patients (Boraud et al. 
[1]). The efliciency of the averaging transformation has been also investigated 
for other brain compartments. In general, the averaged potentials of the com- 
partments reveal oscillatory dynamical behaviour. For example, synchronized 
gamma oscillations were observed in the globus pallidus. The simulation results 
for the globus pallidus could be well reproduced by the reader using the aver- 
aging transformation. Therefore, the representation of these results is omitted 
in the present study. 
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Figure 2: Schematic ultrastructural morphology of the primate's striatum, including the neu- 
rochemical synaptology. The green/orange arrows denote inhibitory/excitatory afferents. The 
dashed arrows are external efferents to the striatal neurons. Hereby, StiA and StiB denote 
the spiny projection neurons which project to Globus Pallidus externa and Globus Pallidus 
interna. St2 denotes the dopaminergic interneurons, St'i the GABAergic interneurons, and 
St4 the acetylcholinergic neurons. 
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Figure 3: Left: The distribution of striatal intcrncurons based on experimental data. Right: 
The oscillatory dynamics of the activity of striatal neurons by a single constant activation 
with a realistic spatial distribution of the neurons (tonically active cholinergic neurons). The 
light nodes denote active neurons. 



4. Discussion 

In conformance with LFP experiments, the multiscale averaging operator (1) 
transforms the synaptic potentials along the network structure of brain region 
such as corpus striatum to local field potentials of the same frequency range. 
This represents a first step from the neuroscicncc at Hodgkin-Huxely dimension 
to the system biological level of consideration. Although this work was initially 
inspired by one of the most intriguing problems of computational neuroscience, 
it appears to be appropriate for several other intriguing applications. The char- 
acterization of the discrete fundamental domains in a complex network in terms 
of n-cells which is one of the dominating concepts of this mathematical ap- 
proach, provides the possibility of the multiscale transformations of information 
along discrete structure that are continuously realized. Beside the advantages of 
this model such as its agreement with electrophysiological experiments, further 
computational investigations are needed to verify its validity for more complex 
systems. 
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Appendix.. 

• Neurotransmitter-NemotrsiBsmitters are the most common class of chem- 
ical messengers in the nervous system; 

• Synapse-Chemical synapses are specialized junctions through which neu- 
rons signal to each other and to non-neuronal cells such as those in muscles. 
Chemical synapses allow neurons to form circuits within the central ner- 
vous system. They are crucial to the biological computations that underlie 
perception and thought; 

• Morphology-The term morphology in biology refers to form, structure and 
configuration of an organism. This includes aspects of the outward appear- 
ance (shape, structure, colour, pattern) as well as the form and structure 
of the internal parts like bones and organs; 

• Ultrastructural morphology-Cell and tissue morphology in the electron- 
microscope level; 

• Intemeuron-An interneuron (also called relay neuron, association neuron 
or local circuit neuron) is a multipolar neuron which connects afferent 
neurons and efferent neurons in neural pathways. Like motor neurons, 
interneuron cell bodies are always located in the central nervous system; 

• Golgi mei/iod-Golgi's method is a nervous tissue staining technique; 
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• Immunohistochemical sfjirfies-Immunohistochcinistry or IHC refers to the 
process of localizing proteins in cells of a tissue section exploiting the 
principle of antibodies binding specifically to antigens in biological tissues; 
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